Synthesizing arbitrary two-photon polarization mixed states 
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Two methods for creating arbitrary two-photon polarization pure states are introduced. Based on 
these, four schemes for creating two-photon polarization mixed states are proposed and analyzed. 
The first two schemes can synthesize completely arbitrary two-qubit mixed states, i.e., control all 15 
free parameters: Scheme I requires several sets of crystals, while Scheme II requires only a single set, 
but relies on decohering the pump beam. Additionally, we describe two further schemes which are 
much easier to implement. Although the total capability of these is still being studied, we show that 
they can synthesize all two-qubit Werner states, maximally entangled mixed states, CoUins-Gisin 
states, and arbitrary Bell-diagonal states. 
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I. INTRODUCTION 

Quantum information processing Q promises great 
power relative to its classical counterpart. Many quan- 
tum information processes require specific pure entan- 
gled states, such as Bell states, to succeed. After in- 
teracting with the environment, however, pure states 
inevitably decohere; decoherence generally causes pure 
entangled states to become mixed and less entangled. 
Quantum error correction Q and entanglement distilla- 
tion/concentration have been developed to help cope 
with a noisy (and hence decohering) environment. On 
the other hand, there are implementations using mixed 
states to investigate quantum computing, e.g., liquid- 
state NMR 01 . The states in this last example are highly 
mixed and have no entanglement. Still, between highly 
entangled pure states and highly separable mixed states 
there exists a vast experimentally unexplored region in 
Hilbert space (more precisely, the space of density ma- 
trices), where states can be simultaneously mixed and 
entangled. The two-qubit system possesses the simplest 
and smallest Hilbert space that permits the existence 
of entanglement. Separate from the specific protocols 
which make use of the states, it is of fundamental in- 
terest to understand the preparation of one of the most 
basic quantum systems. Although there have been many 
attempts _5., 6^, 2., ^ to synthesize two-qubit mixed polar- 
ization states, none has yet been able to create completely 
arbitrary two-qubit mixed states . 

Here we describe several two-photon polarization state 
implementations that should in principle enable prepara- 
tion of arbitrary two-qubit mixed states, including states 
possessing all physical degrees of entanglement and en- 
tropy. The schemes we shall present facilitate state cre- 
ation and allow access to two-qubit Hilbert space and can 



be useful for current and future quantum information 
protocols. We remark that if there exist efficient two- 
qubit entangling gates such as CNOT il| , arbitrary state 
synthesis can be systematically implemented by first gen- 
erating a purification |lfl| | of the mixed state by adding 
ancillas, and then tracing over the ancillas. However, ef- 
ficient photon-polarization CNOT gates do not exist 
so we rely on other degrees of freedom to introduce de- 
coherence, leading to mixed states. 



* Present address; Department of Physics and Optical Engineering, 
Rose-Hulman Institute of Technology, Terre Haute, IN 47803, USA. 



The paper is organized as follows. After a brief back- 
ground discussion in Sec. II Al we describe, in Sec. m 
two schemes to achieve arbitrary pure two-photon po- 
larization states by employing downconversion in a two- 
crystal arrangement. The first one is based on the exis- 
tence of Schmidt decompositions. The second one utilizes 
the coherent superpositions of two downconversion pro- 
cesses embedded in an interferometric setup. In Sec. IIIII 
we describe how to extend these two schemes to real- 
ize arbitrary two-qubit mixed states (Schemes I and II). 
In Sec. IIVI we propose two reduced schemes (HI and 
IV) that provide practical ways to realize several im- 
portant families of states that are currently of inter- 
est, includin g W erner states maximally entangled 
mixed states [T3|. Collins-Gisin states and arbitrary 
Bell-diagonal states. Scheme HI requires only two down- 
conversion crystals, but cannot synthesize all two-qubit 
states. Scheme IV partially extends the set of attainable 
states, but requires four downconversion crystals. Fi- 
nally, in Sec. we summarize the four schemes and men- 
tion possible applications. Readers who do not require 
full details but want an overview of the four schemes and 
synthesizable states can refer to Table |T| The details of 
how to create particular families of states can be found 
around the equations [CHI, OHl), iESIj and (|25J|] describ- 
ing these states. 
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A. Background information 

The entangled photon pairs we consider come from 
frequency-degenerate type-I spontaneous parametric 
downconversion (SPDC) The general state from 

SPDC is a two-mode sc[ueezed state consisting of vac- 
uum and fc-pair states [iq : 

1*) = [vacuum) + e\ijj'^^^) + e^lV-^^^) + • • • , (1) 

where \^^'''>) is a fc-pair state, and e is the relative am- 
plitude (typically of order 10~^) to find a single pair. 
The post-selected 1-pair state IV'^^'') is composed of two 
daughter photons, usually called signal and idler. For 
the present article, we limit our attention to the case 
where the signal and idler photons have approximately 
degenerate central frequencies, half that of the pump. 
(Our schemes apply to nondegenerate case as well.) 
When the downconversion momenta are well coUimated 
or otherwise sharply selected (experimentally by a small 
iris) , one can neglect the momentum dependence of the 
pair state. The post-selected two-photon state can then 
be described by 

\^^^\^)) = {l]c,fc|x,,Xfc)}® f deA,,(e)|| + e,|-e), 



polarization frequency 

(2) 

where w is the pump frequency. \xj^Xk) and \^ +e, ^ —e) 
respectively represent the polarizations and frequencies 
of the two photons, with |xi) = \H), the horizontal 
polarization, and \x2) = the vertical polarization. 
Cjk is the amplitude of the polarization state \xjjXk), 
for single-crystal type-I phase-matching the polarization 
state is unentangled, i.e., cjk = ajbk- ^si(e) is the am- 
plitude for a particular division of energy, so that e indi- 
cates the deviation from half pump frequency. |Asi(e)p is 
peaked at e = with width , and we shall approximate 
it by a gaussian distribution: 

More generally, the pump is not monochromatic, and 
therefore the pair state should be described by 

1^-)= Jdc0Ap{L0)\^/'HL0)), (4) 

where Ap{u}) describes the frequency spread of the pump, 
assumed to be peaked at some frequency loq with half- 
width Suj. For most of the following discussion, we con- 
sider thicknesses of waveplates and crystals that are much 
less than the coherence length Ip (= of the pump, 

and hence we can safely use Eq. (j^J. The coherence 
length of downconversion photons {Isi = c/S^) is usually 
much smaller than Ip, i.e., c/6e <SC c/(5i^, because there are 
many ways to distribute the energy of the pump photon 
between the daughter photons in each pair, resulting in 
a large . 




21PBS I HWP I QWP 

^ NLC Unitary transformation 

FIG. 1: (Color online) Arbitrary pure states via Schmidt 
decomposition. PBS: polarizing beam splitter; HWP: half- 
waveplate; QWP: quarter-waveplate; NLC: nonlinear crys- 
tals. 



II. SCHEMES FOR ARBITRARY 
TWO-PHOTON POLARIZATION PURE STATES 

A. Via Schmidt decomposition 

Using the method of Schmidt decomposition Q], an 
arbitrary two-qubit pure state = a\HH) + b\HV) + 
c\VH) + d\VV) can always be written using only two 
terms:|?A) = a\xO + Plx^^-^), where \x) is orthog- 
onal to Ix^) (1^-^)), and a and /3 satisfy |ap + |/3p = 1. 

Now we describe how to prepare such a state. The 
creation of the entangled state 

cose\HAHB) + e"t' sin^lVAl^s) (5) 

from two-crystal downconversion was proposed in 
Ref. Tsl. Consider two identically cut thin nonlinear 
crystals. Suppose the first crystal's optic axis lies in 
the vertical plane defined by the directions of pump 
beam and the vertical polarization. Assuming type-I 
phase matching, a V^-polarized pump will produce two 
i?-polarized daughter photons. We denote this process 
by 1^^) — *■ \Ha) <8 \Hb)- If the pump is iJ-polarized, no 
downconversion process will take place. Suppose the sec- 
ond crystal is placed at an orientation rotated from the 
first crystal by 90° about the pump direction. An H- 
polarized pump will now produce a pair of 1^-polarized 
photons \H) \Va) \ Vb), whereas no downconversion 
will occur if the pump is F-polarized With the two 
crystals placed in contact with each other, a pump in the 
state cos ^ 1 1^) -I- e*"^ sin 6 \ H) will produce a pair of photons 
in the state 

cose\HAHB)+e"<' sin6l|VAFs), 

where 9 and 4> are tuned using waveplates acting on the 
pump polarization |l5j . {4> can also be tuned with, e.g., 
a variable waveplate acting on just one of the downcon- 
version photons.) 



3 




a\y) + d\H) 
> 



CD 



CD 

Beam Splitter ^ _ 



b\V) + c\H) 



SS 83 



Mirror 



FIG. 2: (Color online) Arbitrary pure states via inteferometry. 

Choosing the local unitary transformations IJ a and Ub 
such that 



{|x>,lx^>}, 

{io,r>}, 



(6a) 
(6b) 



we can achieve the arbitrary two-qubit pure state by 
starting with an entangled state with (cos 0, e**^ sin G) = 
(a, (3) [up to an irrelevant overall phase], followed by the 
corresponding local rotations Ua and U b 

Ua ® Ub{cos9\HaHb) + e"^sin6'|VAV^B)) 

- gIHaHb) + b\HAVB) + c\VaHb) + d\VAVB). (7) 

The two rotations can be obtained in the process of 
Schmidt decomposing \^) 0; see also Appendix IXI for 
an explicit construction of the appropriate C/'s, a, and /3 
given {a, 6, c, d}. 

In practice, any SU (2) rotation such as Ua and Ub on a 
polarization state can be implemented by combinations of 
half- and quarter-waveplates |l9j — preferably zero-order 
waveplates |0, for which the retardance is barely sensi- 
tive to deviation from the central frequency. That is to 
say, the action of waveplates, U, can be assumed to be 
e-independent (at least in the frequency range set by the 
interference filter before detection), i.e., 



J2Uk,\xk)^fdeA{e)\- 



U{\x,)®JdeA{e)\-±e)} 

k 

(8) 

where Ukj are the elements of a unitary matrix that is 
independent of e. We shall assume throughout this paper 
that unitary transformations by waveplates are ideal and 
independent of deviation from the central frequency. 



B. Via inter ferometry 

A second method for creating arbitrary pure states is 
shown in Fig. |21 and is a modification from the setup of 



Ref. |2J. As discussed in Sec. ^1 via SPDC (assuming 
type-I phase matching), a pump in a polarization state 
a\H) + P\V) will generate an entangled photon pair in 
the state (up to some irrelevant phases) a\VV) + /3\HH). 
With a half-waveplate, this entangled state can be fur- 
ther transformed into a\VH) + P\HV). Now, an ar- 
bitrary pure two- photon polarization state a\HH) + 
b\HV) + c\VH) + d\VV) can be regarded as a super- 
position of two (un-normalized) parts: a\HH) + d\VV) 
and b\HV) + c\VH). The first part can be created from 
the (un-normalized) pump state IV'c/) = a\V) + d\H). 
To create the second part, we need the (un-normalized) 
pump state = h\V) +c\H), from which SPDC yields 
the two-photon state b\HH) -\- c\VV). Again, a half- 
waveplate in one arm (or equivalently, passing through 
a quarter-waveplate twice) can transform this state into 
b\HV) + c\VH). By coherently superposing the above 
two processes, as shown in Fig.|21 the fully arbitrary pure 
two-qubit state a\HH) +b\HV) + c\V H) + d\VV) can be 
created. The amplitude of each process, which is deter- 
mined by the relative values of (V'c/IV'c/) and can 
be adjusted by the transmission through the beam split- 
ter. Moreover, coherent superposition can be achieved 
by balancing the two path lengths. Thus, by combining 
a well-balanced interferometer and the process of spon- 
taneous downconversion we can realize arbitrary two- 
photon polarization pure states. In the next section we 
shall describe two schemes capable of producing arbitrary 
two-photon polarization mixed states. 



III. SCHEMES FOR ARBITRARY 
TWO-PHOTON POLARIZATION MIXED 
STATES 

Any two-qubit mixed state can be canonically decom- 
posed as follows U: 



(9) 



where {ji/'i)} are orthonormal eigenstates of p. It is there- 
fore natural to realize p simply by mixing its eigenstates 
with probabilities proportional to their eigenvalues A^. 
As we can synthesize arbitrary pure states from one set 
of crystals, individual synthesis of each is straight- 
forward. 



A. Scheme I: Arbitrary two-qubit mixed states I 

The first mixed-state scheme is shown in Fig. O We 
have four pairs of nonlinear crystals, each generating a 
pure state that, when propagating to the output, arrives 
as \tpi) [l^. There is an attenuator in front of each set of 
crystals (except the first set) such that the pump inten- 
sity li going into the i-th set of crystals is proportional 
to \i (arranged in decreasing order: Ai > A2 > A3 > A4). 
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Schemes 


Synthesizable 
states 


CP 


NLC 


Other 
Optics 


Advantages 


Disadvantages 


T • Pitr R 

1 . rlg.|o 


Arbitrary 
two qubits 


iO 


Q 
O 


OO 


Arbitrary states 


1. Birefringence of crystals causes additional 
rotations and possible decoherence 

2. Requires precise spatial-mode alignment 

3. Narrow opening angles of downconversion 
require long path difference for mixing 

4. Potential loss of downconverted photons 

5. Waveplate imperfection and wedges, esp. 
at early stages, cause beam deviation 


II : Fig.B 


Arbitrary 
two qubits 


15 


2 


48^ 


1. Arbitrary states 

2. Not lossy 

in downconversion 

3. Only two crystals 


1. Requires interferometer stablization 

2. Need to compensate reflection-induced 
transformations from mirrors 

3. Variable beam splitters difficult to tune 

4. Lossy in pump 


III: Fig. in 


Werner 11911. 
Collins-Gisin ^ 
and states 11411 


> 10 


2 


10* 


1. Partially tested 
[6, 33, 34] 

2. Minimal spatial- 
mode matching 

3. Only two crystals 


1. Probably not arbitrary states 

2. No complete theory for more than one 
decoherer per arm 


IV: Fig.E 


States from III, 
Bell-diagonal 
states 13UII 
and states I26II 


> 12 


4 


26 


More states than III 


1. Probably not arbitrary states 

2. Birefringence of crystals causes additional 
rotations and possible decoherence 

3. Requires precise spatial-mode alignment 



TABLE I: Comparison of the four mixed-state schemes. CP stands for controllable parameters (out of 15 in total). The 
nonlinear crystals (NLC) are used in the downconversion process. By "other resources", we include waveplates (where a 
general unitary requires, e.g., 1 half-waveplate and 2 quarter-waveplates, hence counted as 3 elements), mirrors, attenuators, 
prisms, and decoherers, and we assume that the pump is already polarized. The crystal and resource numbers given are sufficient 
to produce all states given in the final column. This resource accounting is intended to indicate the relative complexity of the 
various schemes; however, the numbers listed may be reduced for certain states, or possibly by using clever combinations of 
elements (e.g., reflections which modify polarization). *The resource number listed for Scheme II is several items lower than 
a direct counting from Fig. |3 which was shown for clarity with extra mirrors. "'"The resource number listed for Scheme III is 
counted without pump decoherence and with only a single stage of decoherence, and is thus less than a direct counting from 
Fig. 1^ but is sufficient to synthesize the states listed. 

It is less favorable to attenuate the four downconversion 
pure states to tune the probabilities according to their 
eigenvalues, because direct attenuation of the downcon- 
version photons "would, in general, result in unpaired pho- 
tons, i.e., one of the photons "would be absorbed, but not 
the other 

For a pulsed pump, the mixing is incoherent, as the 
arrival time of the do"wnconversion pair (relative to the 
pump pulse) can, in principle, reveal information on 
"where the pair "was generated. For a CW (continuous 
"wave) pump, one can add a path delay (much greater 
than the pump coherence length [l/^]), to each pair such 
that pair-generation amplitudes at all sets of crystals are 
no longer coherent "with one another. We can thus synthe- 
size p by incoherently mixing its eigenstates "with appro- 
priate "weights. As the do"wnconversion process is much 
more likely to produce one rather than multiple pairs 
[e.g., see Eq. JQ], multiple pairs can be ignored. 



Pump 




I HWP I QWP □ Attenuator 

^ NLC Unitary transformation 

FIG. 3: (Color online) Scheme I employs four sets of nonlin- 
ear crystals. The two-photon state created at the i-th set of 
crystals is chosen such that it is the correct state after 
propagation through the subsequent elements. The necessary 
local unitary transformations at each downconversion location 
can be readily calculated f2^. PBS: polarizing beam splitter; 
HWP: half-waveplate; QWP: quarter-waveplate; NLC: non- 
linear crystals. 



5 



B. Scheme II: Arbitrary two-qubit mixed states II 



The interferometric scheme of Sec. Ill Bi can also be ex- 
tended to create arbitrary mixed states. The full scheme 
is shown in Fig. ^ The coherent superposition method 
of Fig. 12 is used to create each of the four pure states 
\ipi) in the decomposition @ and mix them incoherently 
in proportion to their eigenvalues Ai, as in Scheme I. Ar- 
bitrary weights of mixing can be achieved by controlling 
the transmissions of the beam splitters. In order to mix 
the four parts incoherently, we first use timing informa- 
tion |23| such that the state of the pump is 



i)t. 



(10) 



Here, IV'c/i) a-nd \ijjLi) (both un-normalized) are the two 
parts of the pump state that will, ultimately, yield the 
corresponding pure state |2J|- The factors \i)T 

(i = 1,...,4) encode timing information; there is no 
coherence between paths labelled by distinct values of 
i, i.e., {i\j)T = Sij. For this absence of coherence to 
hold, the path-length difference between any two upper 
(or lower) unmatched paths must be greater than the 
pump coherence length [2^. As long as coherence is 
maintained for the corresponding pairs of states, |V'c/i) 
and \ipLi) (for j = 1, . . . ,4), but the time differences for 
the i's are distinguishable, the output state is the desired 
mixed state, once timing information is traced over, i.e., 
discarded l26l|. Note, however, that there is no differ- 
ence in relative timing between signal and idler photons. 
The timing information is coupled solely to the pump 
photons; because this timing information is traced over 
(ignored), downconversion produces an incoherent mix- 
ture of four two-photon states. Also note that with a 
CW pump, the only possibility to detect any coherence 
in the timing information would be to include similar un- 
balanced interferometers in the downconversion output. 

The difference between the present scheme (II) and 
the previous one (I) is that for Scheme II each of the 
four pure states is created probabilistically in the same 
downconversion crystals, whereas for Scheme I each of 
the four pure states is created in downconversion crystals 
at distinct locations. Both schemes yield arbitrary two- 
qubit mixed states by incoherent temporal mixing. 



IV. RESOURCE-OPTIMIZED SCHEMES FOR 
MIXED STATES 

In this section we describe two reduced schemes (III 
and IV) that provide practical ways to realize several 
important families of states that are currently of inter- 
est. Scheme III, whose feasibility has been demonstrated 
experimentally, emerges as an effort to reduce the num- 
ber of downconversion crystals to two by sacrificing the 
generality of the synthesizable states. Scheme IV further 
extends the set of synthesizable states by employing two 




FIG. 4: (Color online) Arbitrary two-qubit mixed-state syn- 
thesis Scheme II. The transmission probabilities of the var- 
ious beam splitters depend on the desired final state. The 
variable beam splitters immediate preceding the unitary ro- 
tations could also be realized by polarizing beam splitters with 
suitable polarization rotations before and after. 



sets of crystals and the mixing technique introduced in 
Scheme I. 



A. Scheme III: Filling the tangle-entropy plane 

Recall that Scheme I requires the use of, at most, four 
sets of SPDC crystals. Since fewer crystal sets would 
be more economical and likely easier to implement, we 
thus propose a modified scheme, which uses only one set 
of SPDC crystals but relies on "controlled" decoherence. 
Although we do not yet know whether this scheme can 
generate arbitrary two-qubit states, it can synthesize sev- 
eral important families of mixed states, including states 
with all physically allowed values of entanglement (char- 
acterized, e.g., by "tangle") and mixedness (character- 
ized, e.g., by the Hnear entropy 

We use thick birefringent crystals with thickness L as 
"decoherers'.' Their effect on a polarization state of defi- 
nite frequency u is (see, e.g., 0, 28]) 



D(L)|x,)®k)-e 



injLu/c 



\Xj) ^ 



(11) 



where the optic axis is assumed to be along, say, X2, i-e., 
the V direction, and rij is the refractive index for the j- 
th polarization state. The decohering elements entangle 
the polarization and frequency degrees of freedom. In 
the output, only polarizations are detected, so we have 
to trace over the frequency degree of freedom in the joint 
pure state (of polarizations and frequencies). In general. 
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Decoherer 



FIG. 5: (Color online) Scheme III employs decoherence. 
Downconversion photon pairs can be decohered, as well as 
pump photons. The decoherers are thick birefringent crys- 
tals, which separate different polarizations and decrease the 
coherence between them. Also shown is a possible decoher- 
ence on the pump beam: the vertical polarization component 
experiences an adjustable extra delay. 



we are then left with a mixed two-photon polarization 
state. In the present scheme, we can have several deco- 
herers in each arm, along with arbitrary unitary rotations 
between the decoherers (only two are shown in Fig.lSJ. 

In addition to directly decohering the downconversion 
photons, one can also decohere the pump photons before 
downconversion, as indicated in Fig.jS] However, as men- 
tioned previously, the pump typically has a much longer 
coherence length than the downconversion photons do, 
and hence may require much greater relative birefringent 
delays, e.g., unbalanced polarization interferometers, to 
achieve decoherence. 

In the limit we are considering, i.e., Ip 3> |An|Li/2 2> 
Isi (where An = (ny — ^h), and Li and L2 are respec- 
tive thicknesses of the decoherers), decohering the pump 
in addition to the downconversion photons does not pro- 
vide further control beyond simply decohering the down- 
conversion photons. Hence, in the following analysis we 
shall not consider decohering the pump. 

Consider a pure initial polarization state of downcon- 
version pairs (0): IV'^^^) = a\HH) + h\HV) -f c\VH) + 
d\VV) (which is created by the method described in 
Sec. After the decoherers D{Li) and D{L2), one 

in each of the two arms, the state is 

^D{L,)^D{L2)\^^^^). (12) 

Tracing over the frequencies the reduced density ma- 
trix for the polarization state is (with p^, = 1-0) (V'l) 



p = Trc p^ =/ de' 



(13) 



In the limit Li,L2 2> c/(^e|An|), where An is assumed 
to be independent of e, the resulting polarization mixed 
state is (in the {\HH),\HV), \VH), \VV)} basis) 



/ |a|2 fad*\ 



|c|2 

\f*a*d J 



(14) 



where / is a complex function of Li and L2 whose exact 
form depends on ^^^(e). For ^si(e) of gaussian form, as 
in Eq. (jSJ, / is given by 



/ = exp 



An(Li 



c/S, 



-iAn{Li+L2)Lj/2c 



(15) 

Note that |/| < 1, with |/| = 1 for Li = L2. 

The family of two-qubit mixed states described by fTl| 
turns out to be of the form hypothesized by Munro et 
al. in their search for the maximally entangled mixed 
states (MEMS), which define the boundary of physi- 
cally allowed states on the tangle-entropy plane "SQl . 
(The family H14|l also contains other maximally entangled 
mixed states, corresponding to different charaterizations 
of entanglement and entropy "sTl.) Although states in 
this family actually fill the physically allowed region of 
the tangle-entropy plane, this does not mean that the 
family contains all two-qubit mixed states. In fact, the 
family l|14(l has only 4 independent real parameters ex- 
cluding the phase of the off-diagonal element. By includ- 
ing the 6 additional real parameters coming from the two 
arbitrary local SU(2) transformations, we can thus con- 
trol 10 of the 15 real parameters associated with general 
two-qubit mixed states. This assumes a single decoherer 
in each arm. The full capability of Scheme III with an 
arbitrary number of decoherence stages is a difficult the- 
oretical problem that requires further investigation. 

Next we specifically describe how to generate maxi- 
mally entangled mixed states, Werner states, and a par- 
ticular class of mixed states recently discussed by Collins 
and Gisin The maximally entangled mixed states 

found by Munro et al. |l3j are of the form 



Pmems 



Pi(0 = 




for I < r < 1: 



for < r < 





\ 





r\ 

2 




\ 



















1 


U 





k 



(16a) 



(16b) 



Here, an irrelevant phase in the nonzero off-diagonal el- 
ements has been set to zero. For pi{r), we only need to 
generate a pure state of the form 



^,HH)+Vl^\HV) + J-\VV), (17) 



followed by decoherers with thicknesses Li = L2- For 
pii{r), we start with 



\\HH) 



\\HV) + ^\\VV), (18) 



followed by decoherers with thicknesses Li and L2 such 
that |/(Li,L2)| = 3r/2 0. This requires either prior 
knowledge of Asi{e) or the tuning of (Li — L2) so as 
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to obtain the correct reduction factor |/|. Similarly, to 
prepare the Werner states of the form 



pw{r) = r|$+)($+|- 





l-r 
4 









i \ 









(19) 



2 " i±^/ 



[with |<I>^ 
state 



{\HH) + \VV))/V2], we start with the pure 



1- 



\HH) 



l-r 



l-r 



\HV) + \I — \VH) + ^^\VV), 

(20) 

and follow with decoherers with thicknesses Li and L2 
such that |/(Li,L2)| = 2r/(l + r). Analogous procedures 
yield the other forms of the Werner states, i.e., with other 
maximally entangled components. 

Using these methods, several maximally entangled 
mixed states and Werner states have been synthesized 
experimentally, with high fidelities [s^ between the ex- 
perimentally produced states and the theoretical target 
states. For example, the MEMS 









1 
3 






3 







y 



(21) 



was used to investigate entanglement purification proto- 
cols and the separable Werner state 



/i 

1 

6 










1 







(22) 



was used to perform ancilla-assisted process tomography 
without entanglement |8^ . 

Next, we turn to the CoUins-Gisin states, particular 
mixtures of two pure states: 



Pcg(A,0) = Xm{i'e\ + {l-X)\HV){HV\ 



(23) 



A cos^ ( 





Acos0sin6l\ 
(1-A) 




\A cos 6* sin 61 Asin^e* / 

where \ipe) = cos9\HH) +sm9\VV). Collins and Gisin 
reported a Bell-like inequality (which they call 13322) 
that is inequivalent to the usual CHSH-Bell inequal- 
ity [s^ , in that there are states that do not violate CHSH 
but do violate 13322 For example, the family of 

states pcoi^-i^^) exhibit this behavior for certain ranges 
of A and 9 where no violations of CHSH occur. How 



can we create these Collins-Gisin states? In light of the 
above examples of MEMS and Werner states, we see that 
we only need to generate a pure state of the form 

y/Xcos9\HH) + VT^\HV) +VX sine \VV), (24) 

followed by a decoherence with Li = L2. Such states 
have been experimentally realized and used to study var- 
ious tests for entanglement and nonlocality [3^ . 

States described by Eq. (plus those derived from 
them by local unitary transformations) are not the most 
general form that Scheme II can achieve. For exam- 
ple, if, via downconversion, we prepare the pure state 
|*+) = i\HV) + \VH))/y^, apply decoherers of com- 
mon thickness L (^ Igi) in both arms, and then rotate 
each photon polarization by 45°, followed by a second 
set of decoherers with the same thicknesses, we would 
generate a mixed state of the form 



l\ 




1 1 

4 8 

i i 

8 4" 

i 



(25) 



up to some irrelevant phases. This state does not belong 
to the family (|14|) . obtained with only one stage of de- 
coherence, suggesting that using multiple decoherences 
may enable control over more than the 10 independent 
parameters allowed by a single decoherence. Further the- 
oretical investigation is needed to determine the most 
general states obtainable. 



B. Scheme IV: A hybrid technique 

From Scheme I it appears that one needs four sets 
of nonlinear crystals in order to synthesize fully general 
rank-four mixed states, whereas from Scheme III one can 
create rank-four mixed states of the restricted form 1)14(1 
(up to local unitary transformations) with a single set of 
crystals. As we now discuss, by using a hybrid scheme 
one can, with only two sets (or in some cases, even just a 
single set) of crystals, generate a larger class of rank-four 
mixed states. 

The idea is as follows. Suppose that the state p can be 
decomposed into 



(26) 



If the mixed state a can be created by Scheme HI [e.g., 
states in Eq. I14j . we can then mix, with appropriate 
weights, a (as created from a first set of crystals) with 
the pure state \ip){tp\ (from a second set), and thus obtain 
p |3g; see Fig.EK- Although any two-qubit mixed state 
always allows a decomposition into a mixed state plus 
a pure part [3^ . it remains an open question whether 
there always exists a decomposition for which the mixed- 
state part is achievable by Scheme HI. Nevertheless, this 
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(a) 



Pump 












(b) 



FIG. 6: (Color online) Scheme IV is a hybrid technique. 

(a) Mixing a pure state with a mixed state. The local unitary 
transformations immediately after the decoherers are used 
to pre-compensate the effect of local unitary transformations 
used afterward to rotate the pure part, and also to undo any 
effects of passing through the nonlinear crystals (c.f. Fig. 

(b) A reduced setup of the method in (a), using only one 
set of nonlinear crystals, and retro-reflecting the pump back 
through the nonlinear crystals and the first photon pair back 
into the spatial modes of the second pair. Setup (b) is less 
general than (a), as the pure-state part cannot be chosen ar- 
bitrarily. 



hybrid scheme can obviously generate more states than 
Scheme III alone, since adding the pure part adds more 
degrees of freedom. 

We can make a simple reckoning of the number of 
parameters of the achievable density matrices indepen- 
dently controllable. Suppose we restrict the mixed-state 
part a to be produced by Scheme III with only one stage 
of decoherence, i.e., a is of the form H14|l with all param- 
eters real. Recall that a general two-qubit pure state can 
be expressed as 

IV;) =Ua<» Ub{Vx\HV) + Vl - MVH)), (27) 

where Ua and Ub are local unitary transformations, and 
•\/A and \/l — X are Schmidt coefficients. Expressed in 
the basis in which the pure state \tp) is Schmidt decom- 
posed, the mixed state pa + {1 ~ p)\ip) (V'l appears as 



pU{<E>Ul 



+ {1-P) 





^2 








/ fad\ 

62 



\fad 

/O 

A 


\0 



Ua®Ub 







A) 



VA(l-A) 
1 - A 




(28) 



For fixed Ua and Ub, this gives, in general, 6 independent 
parameters {a, 6, c, /, p. A} [noting that d is not indepen- 



dent of {a,6, c}]. Barring some coincidence that, for dif- 
ferent pairs of {a, 6, c, d, /} and {C/4, Ub}, gives the same 
mixed part, we have in total 12 independent parameters, 
after adding 6 from the local unitaries js^. 

One important family of states that this scheme can 
synthesize (and which cannot be generated via Scheme III 
with only one stage of decoherence) are the arbitrary 
Bell-diagonal mixed states Q : 

PB = Ai|$+)($+|+A2|$-)($-|+A3|«'+)(«'+|+A4|*-)(*-|. 

(29) 

Expressed in the {\HH), \HV), \VH), \VV)} basis 



/Ai + A2 



PS =2 











Ai - AsX 



A3 + A4 A3 - A4 
A3-A4 A3 + A4 
VA1-A2 A1+A2/ 



(30) 



Assuming that IA3 — A4I < 1/2 (otherwise |Ai — A2I < 
1/2), pb can also be decomposed as 



A4||*)(«'| 



(31) 



Pb = (1 - IA3 - A4|)pi -I- IA3 
where (1 — IA3 — A4|)pi is 

/A1+A2 Ai-A2\ 

A3+A4-IA3-A4I 

A3+A4-IA3-A4I 

VA1-A2 A1+A2/ 

(32a) 

and the pure-state part is a Bell state 

I*) = ^{\HV)+sgn{X3~\i)\VH)). (32b) 

Here sgn(a;) is the sign function, which gives a factor of 
±1, depending on the sign of x. It is clear that pi be- 
longs to the family of states (|14() . and hence can be syn- 
thesized by Scheme III with one stage of decoherence; on 
the other hand, \'^) is a Bell state, which can be eas- 
ily generated. Furthermore, the weight of pi is not less 
than that of so there is no need to attenuate 

the intensity of the mixed part [s^. Therefore, Scheme 
IV can synthesize any Bell-diagonal state pB- The Bell- 
diagonal states, if entan gled , can be readily distilled via 
the BBPSSW scheme |3ll39| into states with more entan- 
glement. They also have the property that, for a given 
set of eigenvalues, they achieve the maximal violation of 
the CHSH-Bell inequahty 

For certain states this hybrid scheme can also be im- 
plemented via a single set of crystals, by reflecting the 
source and the downconversion pair back through the 
crystals with a mirror; see Fig. ^p. However, in this 
case, the pure-state part cannot be arbitrarily chosen, as 
the local unitary transformations needed to create the 
mixed part a and the pure part 14") are no longer in- 
dependent. The mixed part a is obtained via locally 
rotating ai\HH) + l3i\VV) by IJa ® Ub, followed by 
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a decoherence. As the photons reflect back from the 
mirrors, they experience again the same local unitary 
transformation Ua ®Ub- The local unitary transfor- 
mation at the output port is then chosen to eliminate 
this additional effect (by choosing the inverse of this 
transformation), thereby fixing a to be of the form (|14|l : 
the pure state part is consequently limited to the form 
\^) = IJ-'^ (^IJ-^ {a2\HH) + I32\VV)) . Any further lo- 
cal unitary transformation will rotate a and |^) together 
and cannot change this relative relation. 



V. CONCLUDING REMARKS 

We have described two approaches for synthesizing ar- 
bitary two-qubit pure states. Based on these, we have 
developed four schemes for synthesizing two-qubit pho- 
ton polarization mixed states. Scheme I (Fig.OJ requires 
several sets of downconversion crystals to create arbi- 
trary two-qubit mixed states. It would be desirable to 
experimentally synthesize rank-two mixed states using 
this scheme, in order to give a proof-of-principle demon- 
stration. Scheme II (Fig. 0)) employs temporal mixing to 
achieve decoherence. It offers a second way to realize 
arbitrary two-photon polarization mixed states, but — 
significantly — requires only one set of downconversion 
crystals, at the cost of requiring several rather large, 
phase-stabilized interferometers. Scheme III (Fig.O pro- 
vides control over at least 10 of the independent real 
parameters of two-qubit mixed states, and gives access 
to all physically allowed values of entanglement and en- 
tropy. Furthermore, this scheme has been experimentally 
implemented to synthesize several interesting families of 
mixed states, such as Werner states, maximally entangled 
mixed states, and CoUins-Gisin states 0,113 113 • The 
fourth scheme (Fig. ^ , extends the range of Scheme III 
(providing control over 12 mixed-state parameters). In 
particular, this scheme can be used to produce arbitrary 
Bell-diagonal states, which are of interest, e.g., in en- 
tanglement distillation [J and maximal violations of 
Bell inequalities 'i^. Although the full capabilities of 
Schemes III and IV are not yet entirely clear, our analysis 
shows that these two schemes provide practical methods 
for creating quite general mixed states, many of which 
were previously not accessible experimentally. The four 
mixed-state schemes are summarized and compared in 
Table including the respective resources, advantages 
and disadvantages for implementation. 

Once one has well-controlled arbitrary two-qubit 
sources, they will be usable for many quantum infor- 
mation processing applications, such as testing meth- 
ods of entanglement distillation 0, lEJsMEJl' investi- 
gating quantum process tomography |33l l43 | , character- 
izing quantum gates [431. testing violations of Bell-type 
inequalities d, 0, l35l l40| by mixed states (including a 
relevant two-qubit Bell inequality [T3 | that is inequiva- 
lent to the CHSH inequality), and exploring the vast, 
previously inaccessible territory of Hilbert space. 
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APPENDIX A: CREATING ARBITRARY PURE 
STATES 

In this appendix we explain how to create an arbi- 
trary two-qubit pure state, characterized by {a, 6, c, d\ of 
Eq. ^ . This amounts to establishing adequate local uni- 
tary transformations {Ua and Ub) and post-SPDC pure 
states Eq. of the form that SPDC naturally yields. 
For convenience, we exchange the coefficients cos for a 
and e*"^ sin 6 for /3. 

To find the appropriate settings and local unitary 
transformations, we need to solve 

UA®UB[a\HH)+(}\VV)) = 
a\HH)^h\HV) + c\VH)+d\VV), (Al) 

for {Ua, Ub, ot, /?}, given {a, 6, c, d\ that are properly nor- 
malized. This equation can be solved either by Schmidt 
decomposition or by direct algebraic manipulation. How- 
ever, the solution is not unique. 

When ad—bc— 0, the state to synthesize is a product 
state, which can be created from an initial state \HH) fol- 
lowed by independent local rotations (see, e.g., Ref. (igj). 
For 

ad-bc^O, \ad-bc\^ 1/2, (A2) 

i.e., the case of non-maximally entangled pure states, one 
possible solution is 

a = \Jl- V'l~4|ad-6c|V%/2, (A3a) 
13= {ad- be) /a, (A3b) 

= { ""l ^i) ^ Ub=( , (A3c) 

where 

ui EE Izil/Vkip + Nap, vi = Z3ul/zl, (A4a) 

U2 = Zy{ziZ3 + Z2ZA), V2 EE Z2 / {ziZ^ + Z2Z4), (A4b) 

with Zl, Z2, Z3, Z4 defined via 

Zl = ( aa*-dy)/{\af~\/3f), (A5a) 

Z2 = { d*a-a/3*)/(|a|2-|/3|2), (A5b) 

^3 ^ (_c*a-6/3*)/(H2-|/3|2), (A5c) 

^4 = ^_ba*-e*f3)/{\a\^-\(3\'). (A5d) 
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When \ad — bc\ = 1/2 there are three possible cases: 

(i) 6 = c = and \a\ = \d\ = 1/^2; 

(ii) a = d = and |6| = |c| = l/\/2; 
(hi) a,b,c,d^O, \a\ = Ml, |5| = |c|. 

Case (i) is aheady the form we seek. In case (ii), an 
exchange H ^ V (e.g., by a half-waveplate) will do. In 
case (iii), one possible solution is 



e 

(a,/3) 



(e^M)/V2, 



(A6a) 
(A6b) 



In fact, this last case includes the previous two cases, if 
one interprets the phase in Eq. HA6a|l to be the appro- 
priate ratio of the nonzero coefficients. We remark that 
the solution presented above is not unique, and that ex- 
perimentally one would implement the one that is most 
convenient. 



Ua = ^ 



Ub = 



V2 \-e-'-^ 1 I 

' d* -c d* +(^ 
-d~c* d-c*. 



(A6c) 
(A6d) 
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